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We present a lifetime measurements of the 6s level of rubidium. We use a time-correlated single- 
photon counting technique on two different samples of rubidium atoms. A vapor cell with variable 
rubidium density and a sample of atoms confined and cooled in a magneto-optical trap. The 5P\/2 
level serves as the resonant intermediate step for the two step excitation to the 6s level. We detect 
the decay of the 6s level through the cascade fluorescence of the 5P3/2 level at 780 nm. The two 
samples have different systematic effects, but we obtain consistent results that averaged give a 
lifetime of 45.57 ±0.17 ns. 

PACS numbers: 32.70.Cs, 32.80.Pj, 32.10.Dk 



I. INTRODUCTION 

This paper presents our measurements of the lifetime 
of the 6s level in rubidium using time-correlated single- 
photon counting techniques in two different atomic envi- 
ronments: An atomic vapor cell, and a magneto-optical 
trap (MOT). This work complements and enhances our 
program of francium spectroscopy and weak interaction 
physics. We have previously measured lifetimes in fran- 
cium and rubidium only in a 

use the same apparatus except for the source of atoms so 
we carefully address the different systematic effects that 
are unique to each one of them. The present measure- 
ments are necessary to understand systematic effects on 
the measurement of the equivalent level in francium, the 
8s level. 

Measurements of excited state atomic lifetimes in the 
low-lying states of the s manifold enhance our under- 
standing of the wave functions and the importance of cor- 
relation corrections in the theoretical calculations. Rela- 
tivistic corrections in rubidium are smaller than in other 
heavier alkali such as cesium or francium. This particular 
level (6s) is of primary importance to the optical PNC 
measurements that look for the parity forbidden dipolc 
transition between the ground state and the first excited 
state with the same orbital angular momentum.^, 
The comparison of measurements with theoretical pre- 
dictions test the quality of the computed wave functions. 
The calculation of the wave functions have now reached 
new levels of sophistication [Tol ITl| based on Many Body 
Perturbation Theory (MBPT). Those calculations are 
particularly important in the interpretation of precision 
tests of discrete symmetries in atoms such as Parit y N on- 
Conservation (PNC) @ and Time Reversal (TR).[T^| 

The paper is structured as follows: We show the rela- 
tionship between lifetimes and atomic structure in sec- 
tion [n] Sec. IIIII presents the experimental method for 
the lifetime measurement with details on the two sources 
of atoms. Sec. IIVI has the data analysis for the life- 
time measurement. Sec. Ivl summarizes the results and 



compares them with previous results and theoretical pre- 
dictions. Sec. IVII contains the conclusions in the context 
of similar measurements in rubidium and francium. 



II. LIFETIMES AND ATOMIC STRUCTURE 

The lifetime of an exited atomic state depends on the 
initial and final state wave functions and the interaction 
that connects them. Since the electromagnetic interac- 
tion in atomic physics is well understood, radiative life- 
times give information on atomic structure. 

The lifetime r of an excited state is determined by its 
individual decay rates, 1/tj, to other states by 



T ^ Ti 

The matrix element associated with a partial decay rate 
(i) between two states connected with an allowed dipole 
transition in free space is: 



i_4or3 |(J||r||J')| 2 

n ~ 3 c 2 2 J' + 1 ' ( ' 

where lo is the transition energy divided by h, c is the 
speed of light, a is the fine-structure constant, J' and 
J are respectively, the initial and final state angular 
momenta, Tj is the excited state partial lifetime, and 
I (J 1 1 r 1 1 J') I is the reduced matrix element. 

A calculation of the reduced matrix element requires 
the precise knowledge of the electronic wave functions in- 
volved. The presence of the radial operator makes the 
matrix element more sensitive to contributions of the 
wave functions at large distances from the nucleus. 
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FIG. 1: Energy levels of Rb. The figure shows the two step 
excitation (solid thick lines), the fluorescence detection used 
in the lifetime measurement (dashed line) and the undetected 
fluorescence (dotted line). The trap and repumper beams 
(solid thin lines) are only used for the measurement in the 
MOT. 



III. EXPERIMENTAL METHOD 

We use time-correlated single-photon counting |l3j to 
obtain the lifetime of the 6s level in rubidium in a vapor 
cell and in a MOT. This two step transition method is 
well established both for vapor cells [l4| and MOTs 0, @ 
in alkalis. 

Figure ^ shows the energy levels of 85 Rb (J = 5/2) 
relevant to the lifetime measurement. We excite the 6s 
level in a two-photon resonant process. The first step 
laser at 795 nm populates the bP\/ 2 F = 3 level and the 
second step laser at 1324 nm excites the 5Pi/ 2 — > 6S1/2 
transition. 

After we turn off the excitation lasers, the atom re- 
turns back to the ground level using two different decay 
channels. First, by emitting a 1324 nm photon it decays 
back to the 5Pi/2 state and fluoresces 795 nm light to 
return to the 5s ground level. This decay path is not 
used for the 6s level lifetime measurement. The second 
possible decay channel is the 6s — ■> 5P3/2 transition fol- 
lowed by the decay to the 5s ground level. The 1367 
nm fluorescence from the first step of this decay is unob- 
served, but we detect 780 nm light from the second part 
of the deca/y. Since the lifetime of the 5P3/2 state is well 
known. |2l Il5l ll6L Il7j it is possible to extract the 6s level 
lifetime indirectly from the fluorescence distribution of 
the second step of the cascading decay. 

We use a high efficiency MOT to capture a sample of 
rubidium atoms at a temperature lower than 300 /LtK. 
The MOT is in line with the Stony Brook Superconduct- 
ing LIN AC and is optimized to trap francium. |Tsj We 
load the MOT from a rubidium vapor produced by a dis- 
penser in a glass cell coated with a dry film to reduce 
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FIG. 2: Experimental setup. The trap and repumper beams 
are not used for the vapor cell measurement. The dotted 
(dashed) boxes are not used for the MOT (vapor) experiment. 
AM EOM=amplitud modulated electro-optic modulator, AM 
AOM=amplitude modulated acousto-optic modulator 



sticking to the walls. The MOT consists of three pairs of 
retro-reflected beams, each with 15 mW/cm 2 intensity, 
3 cm diameter (1/e intensity) and red detuned 19 MHz 
from the atomic resonance. A pair of coils generates a 
magnetic field gradient of 9 G/cm. The trap contains 10 4 
atoms with a diameter of 0.5 mm and a typical lifetime 
between 5 and 10 s. 

We perform the lifetime measurement also in a vapor 
cell. The cell is under vacuum (approximately 10~ 5 Pa) 
and a dispenser provides the rubidium for the measure- 
ment. The cell is uncoated and we use no buffer gas. We 
keep the rubidium density low to avoid collisional quench- 
ing and radiation trapping effects. The typical mean free 
path of the atoms in the cell is more than 30 m. We can 
apply a uniform magnetic field to the cell in contrast to 
the permanent magnetic field gradient in the MOT. 



A. Lasers System 

Figure 21 shows the experimental setup. We use a Co- 
herent 899-21 titanium-sapphire (Ti:sapph) laser at 780 
nm between the 5s F = 3 and the 5P3/2 F — 4 levels 
to trap and cool the atoms. A Coherent 899-21 Thsapph 
laser at 795 nm between the 5s F — 2 and the 5Pi/ 2 
F = 3 repumps any atom that falls out of the cycling 
transition. The trap and repumper beams are not nec- 
essary for the vapor cell measurement. We excite the 
atoms to the 6s level via a two step transition. A Co- 
herent 899-01 (or a Coherent 899:21 for the vapor cell) 
at 795 nm between the 5s F = 3 and the 5Pi/2 F = 3 
levels makes the first step to the 5P±/2 level and a EOSI 
2010 diode laser at 1.324 fim completes the transition to 
the 6s level. 

A Burleigh WA-1500 wavemeter monitors the wave- 
length of all the lasers to about ±0.001 cm -1 . We lock 
the t rap , repumper and first step laser with a transfer 
lock. [l3 We use a different method to lock the second 
step laser on the vapor cell and on the MOT measure- 
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merits as described below. Figure shows both locking 
methods for the second step laser. 

Second step laser lock for the MOT. We transfer 
the stability of the Helium-Neon laser (Melles-Griot He- 
Ne 05-STP-901) used on the lock of the other lasers 
to lock a Michelson interferometer with a 200 MHz free 
spectral range. The 1324 nm light propagates collinearly 
with the He-Ne light inside the interferometer. We detect 
the output of the locked interferometer at 1324 nm and 
use that signal to lock the second step laser to a side 
of the fringe, which allows us to set it on the atomic 
resonance by an offset adjustment. We achieve a drift 
smaller than 5 MHz per hour. 

Second step laser lock for the vapor cell. We use 
two photon modulation transfer spectroscopy to lock the 
second step laser. We split the first step laser at 795 nm 
into two beams and send both of them through a sepa- 
rate rubidium glass cell. The room temperature vapor 
cell is uncoated and the beam propagates approximately 
7 cm inside the cell, which leads to a typical absorption 
of about 15 %. We collinearly apply the 1.324 fim second 
step laser to one of the beams to excite the second tran- 
sition. If the laser is resonant to the two step transition 
it results in less absorption of the 795 nm light. Two 
beams and balanced detection allow common mode noise 
rejection to increase the signal to noise ratio. We use the 
internal piezoelectric crystal of the 1.324 /im diode laser 
to dither its frequency at 1.3 kHz and perform lock-in 
detection of the 795 nm light transmission through the 
cell. We feedback the resulting error signal to the second 
step laser to lock it with comparable performance to that 
of the Michelson lock. 

We send and retro-reflect 1.5 mW of the first step laser 
power to the trap region. We focus it to an spot size 
approximately equal to the trap size. The second step 
laser travels at an small angle with respect to the first 
step laser and we also focus it to match the trap size 
with a power of 1 mW. For the vapor cell measurement 
the excitation beams propagate collinearly inside of the 
glass cell. We send 100 [iW of power at 795 nm and 0.4 
mW of power at 1.324 /mi and focus it to a spot size of 
about 0.6 mm. 
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FIG. 3: Block diagram for the 1324 nm second step laser 
locking system for the measurement: a) with the MOT and 
b) with the vapor cell. 
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FIG. 4: Timing diagram for the 6s level excitation and decay 
cycle. The cycle repetition rate is 100 kHz. There is no trap 
nor repumper for the vapor cell measurement. 



B. Timing 

Figurc0]displays the timing sequence for the excitation 
and decay cycle. We apply the two step excitation 500 ns 
after the trap laser turnoff for a duration of 100 ns. We 
turn off the repumper at the same time as the two step 
lasers to look for the fluorescence decay of the atoms. The 
counting electronics are sensitive for 500 ns to record the 
excitation and decay signal. The cycle repetition rate is 
100 kHz. The trap and repumper pulses are not necessary 
for the vapor cell measurement. 

Modulation for the MOT. We turn the trap laser on 
and off with an electro optic modulator (EOM) (Gsanger 
LM0202) and an acousto-optic modulator (AOM) (Crys- 



tal Technology 3200-144) as shown in Fig. H The EOM 
gives a fast turnoff and the AOM improves the long term 
on/off ratio. The combination of the two gives an extinc- 
tion ratio of better than 1600:1 after 500 ns. We mod- 
ulate the repumper and first step laser with an AOM 
(Crystal Technology 3200-144) with an on/off ratio of 
109:1 and 26:1 respectively after 30 ns. We couple the 
1324 nm light into a single mode optical fiber. We use a 
combination of a 10 Gbits/s lithium niobate electro-optic 
fiber modulator (Lucent Technologies 2623N) followed by 
an fiber amplifier (InPhenix IPSAD1301) and a second 
identical fiber modulator to turn the laser on and off. 
The combination of the three elements gives an on/off 
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ratio better than 1000:1 after 20 ns. 

Modulation for the vapor cell. For the vapor cell 
measurement we do not need the trap and repumpcr 
beams and we install the EOM on the first step laser 
path instead. In this case we achieve an extinction ratio 
of better than 100:1 in 30 ns. We reduce any radio fre- 
quency (RF) emission generated by the fast turnoff of the 
EOM placing it in a separate room inside a metal cage. 
We use only the last fiber modulator for the 1324 nm 
laser and we achieve an on/off ratio of 110:1 after 20 ns. 
Figure |21 shows an schematic diagram of the modulation 
of all the lasers. 



C. Imaging system 

A 1:1 imaging system (f/3.9) collects the MOT fluo- 
rescence photons onto a charge coupled device (CCD) 
camera (Roper Scientific, MicroMax 1300YHS-DIF). We 
monitor the trap with the use of an interference filter at 
780 nm in front of the camera. A beam-splitter in the 
imaging system sends 50 % of the light onto a photomul- 
tiplier tube (PMT) (Hamamatsu R636). An interference 
filter at 780 nm in front of the PMT reduces the back- 
ground light other than fluorescence from the cascade 
through the 5P3/2 level decay back to the ground state 
55 1/2 . 

The 1:1 (f/2.4) imaging system for the vapor cell is 
perpendicular to the excitation beams and collects the 
fluorescence from the decaying atoms to a PMT (Hama- 
matsu R636) with a 780 nm interference filter in front of 
it. 



D. Detection electronics 

Figure[S]shows the electronic processing of the detected 
fluorescence. When the PMT detects a photon, it gen- 
erates a variable amplitude electronic pulse. We first 
amplify the signal with an Ortec AN106/N amplifier, we 
send it to a linear gate (EG&G LG101/N, not used for 
the vapor cell measurement) and then to an Ortec 934 
constant fraction discriminator (CFD). The output of the 
CFD goes to a gated time-to-amplitude converter (TAC) 
(Ortec 467). The gating is done directly on the TAC for 
the MOT measurement and with an external gate (Ortec 
LG101/N) for the vapor cell measurement. The photon 
pulse starts the TAC and we stop it with an electronically 
generated pulse that has a fixed time delay with respect 
to the excitation laser pulse. Starting the TAC with a flu- 
orescence photon eliminates the accumulation of counts 
from cycles with no detected photons. We use a multi- 
channel analyzer (MCA) (EG&G Trump-8k) to produce 
a histogram of the events showing directly the exponen- 
tial decay. A Berkeley Nucleonics Corporation BNC 8010 
pulse generator is the master clock with slaved pulse and 
signal generators to provide the timing sequence for the 
measurement. 
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FIG. 5: Block diagram of the processing electronics. 
The dashed gate is not used for the vapor cell experi- 
ment. PMT=photomultiplier tube, CFD=constant frac- 
tion discriminator, TAC=time to amplitude converter, 
MCA=multichannel analyzer. 



IV. DATA ANALYSIS 

We extract the atomic radiative lifetime from the mea- 
sured fluorescence decay. We take sets of data for about 
600 s, that are individually processed, added together, 
and fitted. The 6s level decays through the 5p states to 
the 5s level as shown in Fig. and we detect the indirect 
decay through the 5P 3 / 2 state. The fitting function is a 
sum of two exponentials, a background and an sloping 
background coming from the trap laser turnoff 



A 6s exp 



T6s 



A 5p exp ( -— ) (3) 



-A R + Ast, 



where t§ p is the known lifetime of the 5P 3 / 2 state and tq s 
the lifetime we want to extract, Aq s and A$ p are the coef- 
ficients of the exponentials while Ag and Ag characterize 
the background. Ag = for the vapor cell measurement 
since there is no trap laser. Equation 0] reduces to a sin- 
gle exponential if one waits long enough after the turnoff 
to take data. [13 We instead maximize the statistics by 
including both exponentials in our analysis. 

The first part in the data analysis requires a pile-up 
correction, that accounts for the preferential counting of 
early events. |13| If iVj is the number of counts in the 
MCA channel i, and Ne is the total number of excita- 
tion cycles, then N- is the corrected number of counts in 
channel i given by 



N' 



1 - — y 

1 N E ^ 



=1 Iy j 



(4) 



As low count rates keep this correction small, we collect 
data with a small number of fluorescence photons. We 
typically count one photon every 100 cycles, which cor- 
responds to a correction in the fitted lifetime by 0.1 %. 

We perform a nonlinear least square fit using 
Microcal™ Origin™ Version 5.0 to find the fitting pa- 
rameters that produce the smallest \ 2 ■ We have tested 
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FIG. 6: Decay curve of the 6s level through the 5P3/2 state 
in the vapor cell with fit and residuals. The upper plot shows 
a) the raw arrival time histogram data for the vapor cell mea- 
surement and b) the data after the pile-up correction and 
substraction of background and the 5P3/2 decay. The lower 
plot shows the normalized residuals. Xi5=0.98 



the software independently to validate and understand 
its results. Figure [S] shows an example of the fitting pro- 
cedure. We start the fit 20 ns (25 ns for the vapor cell 
measurement) after both excitation lasers are off. The 
fitting function describes the data well, and the reduced 
xl of this particular decay is 0.98. The average xl f° r the 
MOT data is 1.05 ± 0.06 and for the vapor cell 1.09 ± 
0.08. A discrete Fourier transform of the residuals shows 
no structure. 



The MOT and vapor cell measurements have some 
characteristics that are intrinsically different between 
them. The MOT cools the atoms to a very low tempera- 
ture, while the atoms in the vapor cell have temperatures 
corresponding to Doppler shifts significantly larger than 
the natural linewidth of the transitions of interest. The 
MOT has magnetic field gradient always present which 
can be removed for the vapor cell. These differences gen- 
erate systematic effects that are intrinsic to each system. 
We classify the systematic effects into two groups, com- 
mon effects to both the MOT and the vapor cell mea- 
surement and distinct effects that appear in a different 
way in both systems. 



A. Common systematic effects 

1. Truncation uncertainty 

The truncation uncertainty takes into account the vari- 
ation of the fitted lifetime depending on the start and 
end points of the fit. We find no statistically significant 
change in the obtained lifetime when we change both the 
initial or final point of the fit. 

2. Time calibration 

We calibrate the scale on the MCA by sending a se- 
ries of artificial start and stop pulses with known delay 
directly to the TAC either from one of the DG535 pulse 
generators or from two DG535 pulse generators triggered 
by a frequency synthesizer (Rockland 5600). The linear 
fit to the resulting data contributes a ± 0.01 % to the life- 
time uncertainty and it is consistent using either source 
of calibration pulses. 

3. TAC and MCA response nonuniformity 

A nonuniformity in the TAC or MCA results in differ- 
ent bins of the MCA having different response sensitiv- 
ities that alters the obtained lifetime. We measure the 
height uniformity of the TAC/MCA system by trigger- 
ing the PMT with random photons from a white light 
source. The nonuniformity of the TAC/MCA contribute 
a 0.11 % (0.02 % for the vapor cell measurement) to 
the uncertainty of the lifetime. The need for additional 
electronic gates in the MOT measurement increases the 
nonuniformity in the electronics response. The response 
changes slightly between days. It shows up in the data 
as an additional noise contribution that is responsible for 
the deviation of the xl from one. The effect is already 
included in this uncertainty contribution. 

4- Imperfect laser turnoff 

An imperfect turnoff of the excitation lasers influences 
the observed decay signal. If the second step laser is not 
completely turned off during the decay, atoms that decay 
to the 5Pi/ 2 state can be reexcited to the 6s level. We 
test the turnoff by comparing the lifetime obtained when 
we leave the first step laser continuously on or off during 
the decay. We establish a limit of 0.07 % for the effect 
of the imperfect turnoff of the second step laser in the 
MOT measurement. 

We use two modulators for the second step laser on the 
MOT measurement to achieve a large on/off ratio. The 
vapor cell measurement uses only one modulator and we 
apply a correction to the obtained lifetime. 

The appendix presents the model used to calculate the 
effect of an imperfect turnoff on the lifetime. We solve 
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the rate equations and compare the result obtained with 
a perfect turnoff to the experimental one. We find a cor- 
rection to the lifetime of - 0.11 ns and a sensitivity anal- 
ysis of the parameters of the model gives an uncertainty 
contribution of ± 0.1 %. 



5. Initial state conditions 

The lifetime should be independent of the populations 
of the different levels at the beginning of the decay. We 
can change the initial state conditions by changing the 
power of the first or second step laser pulses, or their 
duration. We calculate the linear correlation coefficient 
and its integral probability to study the correlation be- 
tween the obtained lifetime and the external variable. 
The larger the integral probability is, the less likely it is 
that a correlation exists 

We vary the power of the first and second step lasers 
by more than an order of magnitude. In the vapor cell 
measurement we also changed the excitation pulse dura- 
tion between 50 ns and 200 ns. In all the above cases we 
found results consistent with no correlation. 



B. Distinct systematic effects 

1. Radiation trapping 

The fluorescence photons from atoms decaying from 
the 5P 3 / 2 state to the 5s ground level can excite other 
atoms that are in the ground state and take them to 
the 5P 3 /2 level, the so called radiation trapping. [22I |23| 
Traps are very convenient systems to make measurements 
because the atoms are collected into an small region of 
space. The effect of radiation trapping is small in a low 
density trap because the distance d corresponds to the 
diameter of the trap which is rather small. Vapor cells 
on the other hand may have a large density, but Doppler 
shifts reduce the effective number of atoms that could 
reabsorb the fluorescence. 

The probability of a photon being absorbed by an atom 
can be deduced from Beer's law 



1(d) = Iq exp(— ad), 



(5) 



where Jo is the initial intensity and 1(d) is the intensity 
after traversing the distance d and a — an represents the 
absorption coefficient, with n the density of the atomic 
vapor and a = 3A 2 //2-7r the absorption cross section. 
/ = 0.69 is the oscillator strength of the corresponding 
transition. For small absorption the probability P of a 
photon being absorbed by the atoms is P — an d. 

The rate equation for the 5P 3 /2 level in the presence 
of radiation trapping after the lasers turnoff is given to 
first order by 



dN ; 
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(6) 
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FIG. 7: Temperature dependence of the measured 6s level 
lifetime together with a fit. 



with N 5p and N es the number of atoms in the 5P 3 / 2 and 
6s levels and B$ p the branching ratio from the 6s level 
to the 5P 3 / 2 - In this approximation the change in the 
lifetime of the 5P3/2 level is directly proportional to P. 
As our fitting process assumes a fixed 5P 3 / 2 state lifetime, 
the change due to radiation trapping propagates linearly 
to the measured 6s level lifetime. 

We vary the density of rubidium atoms to look for vari- 
ations in the lifetime. We change the density by chang- 
ing the number of atoms in the MOT or by compressing 
the trap with an increase the magnetic field gradient. 
We establish a limit of 0.01 % for the effect of radiation 
trapping on the lifetime. 

A rubidium dispenser (SAES) provides the atoms in 
the vapor cell measurement. A zirconium alloy matrix 
releases rubidium atoms when it is heated by a DC cur- 
rent. The desorption follows an Arrhenius process that 
requires a certain activation energy E a . The atomic den- 
sity in the cell follows the relation 



n(T) oc exp 



k B T 
E a 



(7) 



where ks is the Boltzmann constant and T the tempera- 
ture in K. Measurements of the pressure on the cell con- 
firm this relationship. 

We compare the lifetime obtained at different dispenser 
temperatures and extrapolate the observed dependency 
to zero temperature to obtain the correction to the life- 
time. The fitting function for the temperature depen- 
dence of the observed lifetime is 



t(T) = To +Aexp(T/B), 



(8) 



where A, B and the corrected lifetime To are fitting con- 
stants. Figure shows the measured data and the ob- 
tained fit. The correction of the 6s level lifetime due to 
radiation trapping amounts —0.28 ns and generates an 
uncertainty contribution of 0.38 % becoming the domi- 
nant uncertainty contribution in the vapor cell measure- 
ment. 



7 



2. Magnetic field 

The presence of a magnetic field can influence the mea- 
sured lifetime through quantum beats between Zeeman 
sublevels. In ideal conditions the atoms in a MOT sit 
at the zero of the magnetic field gradient but in reality 
an imbalance of the power of the retro-reflected beams 
may displace the trap from this ideal position. The use 
of a cascade decay reduces the appearance of quantum 
beats. [Hill 

We quantify the magnetic sensitivity of the lifetime 
measurement in two ways. First we perform the mea- 
surement at different magnetic field gradients. We do not 
observe a correlation with the magnetic field gradient and 
we establish a limit on the effect of the magnetic field on 
the lifetime of 0.17 %. This limit is consistent with mea- 
surements done in the vapor cell. Alternatively we keep 
the magnetic field environment fixed and we displace the 
trap by creating an imbalance on the laser beams that 
interact with it. We can imbalance the trap beams by 
inserting a piece of glass in front of one of the retro- 
reflection mirrors, and repeating the same for all three 
axis. Changing the power of the first step laser beam 
also displaces the trap. The change in the obtained life- 
time in all cases is consistent with statistical fluctuations. 

The measurement in the vapor cell does not require 
a magnetic field gradient. We generate an homoge- 
neous magnetic field with two coils arranged close to a 
Helmholtz configuration with its axes in the same direc- 
tion as the excitation laser beams. The magnetic field 
homogeneity in the excitation region is ± 0.35 %. 

Figure |S| shows the measured 6s level lifetime in the 
vapor cell as a function of the external uniform mag- 
netic field. We can control the magnetic field to ±0.2 
Gauss due to the earth magnetic field and other mag- 
netic sources. The correlation coefficient for the three 
data points is consistent with no correlation. The dotted 
line shows the 68% confidence band (la) on the linear 
fit for a t distribution with one degree of freedom. The 
dashed line gives the error on the average assuming no 
linear dependence. The figure shows that both regions 
are consistent with each other. We limit the contribu- 
tion to the uncertainty due to magnetic fields to ± 0.07 
% using the linear fit to the data. 
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FIG. 8: Magnetic field dependence of measured lifetime in the 
vapor cell. The plot contains the la region for the linear fit to 
the data (dotted line) and the error on the average assuming 
no linear dependence (dashed line) . 



4- Time of flight 

A moving atom can escape the imaging area before 
it decays. The effect is more evident in the vapor cell 
where the atoms have larger speeds. The MOT (vapor 
cell) imaging area has a 1 mm (5.5 mm) diameter aper- 
ture. It takes an atom on average 10 5 (200) lifetimes to 
leave the MOT (vapor cell) imaging area. A Monte-Carlo 
simulation shows no statistically significant change in the 
lifetime for the vapor cell measurement. 



5. Background slope 

The detected light at 780 nm comes from the cascade 
decay from the 6s level, but also from scattered light from 
the trap laser. We turn off the trap laser 500 ns before the 
two step excitation, but the remanent of the trap laser 
light appear in the signal as an slope in the background. 
Fits to files with and without decaying atoms give con- 
sistent values for the background slope. We compare the 
lifetime obtained when we leave the background slope as 
a free parameter or when we fix it to the value obtained 
without two step excitation and obtain an uncertainty 
contribution of ± 0.22 %. This contribution does not ap- 
pear in the vapor cell measurement since there is no trap 
laser. 



3. Collisional quenching 

Inelastic collisions can modify the lifetime of an excited 
state. The velocity of atoms in the MOT (vapor cell) is 
smaller than 0.1 m/s (250 m/s). The mean free path of 
the atoms in the MOT (vapor cell) is about 10 4 m (10 
m) the effect of collisions is negligible for both samples. 



6. PMT response 

We look for saturation effects in the PMT due to the 
trap light by comparing the PMT response to that of a 
fast photodiode which is not subject to saturation. Wc 
replaced the two step excitation by a light pulse at 780 
nm. Comparison of the detected pulse in the PMT and 
in the fast photodiode gives a maximum contribution to 
the uncertainty of ± 0.21 %. This effect is not present in 
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the vapor cell measurement. 



V. RESULTS AND COMPARISON WITH 
THEORY 

The lifetime of the 6s level depends on the value of the 
lifetime of the 5P3/2 state. The uncertainty of the 5P3/2 
state lifetime propagates to the uncertainty of the 6s level 
lifetime. Using Bayesian statistics the contribution to the 
uncertainty (<tb 6> ) due to the 5P3/2 level is given by|| 



Correction % error 



Time calibration 




±0.01 


Bayesian error 




±0.11 


TAC/MCA nonuniformity 




±0.02 


Imperfect laser turnoff 


-0.11 ns 


±0.10 


Radiation trapping 


-0.28 ns 


±0.38 


Magnetic Field 




±0.07 


Statistical error 




±0.24 


Total 




±0.48 



dT 6s (T') 



(9) 



We assume different values r' for the lifetime of r 5 p 3/2 and 
include them in the fitting function (Eq. to obtain the 
dependence of t 6s (t'). 

The 5P3/2 level lifetime has been measured previously 
by several groups (see Table HJ. We take the weighted 
average of all the measurements. We add the systemat- 
ical contributions of the respective errors in quadrature 
and calculate the total statistical contribution according 
to (Tg tat — EV°f) 1 where Oi is the statistical un- 
certainty of a single result. The 5P 3 / 2 state lifetime of 
26.23(9) ns gives a Bayesian error for the 6s level lifetime 
of 0.17 % (0.11 % for the vapor cell). 



Author 



lifetime [nsl 



Simsarian et al. [2] 26.20(9) 

Volz et al. [15] 26.24(4) 

Gutterres et al. [16] 26.25(8) 

Heinzen et al. [17J 26.23(6) 

TABLE I: Experimental results for the rubidium 5P3/2 life- 
time 



Tables ITTI and IIIII contain the error budget and the cor- 
rections for the 6s level lifetime measurement in the va- 
por cell and in the MOT. The error is dominated by the 
statistical uncertainty in the MOT and by radiation trap- 
ping in the vapor cell. The vapor cell corrections change 
the result by -0.9 %. We obtain a lifetime of 45.64 ± 
0.22 ns in the vapor cell and 45.48 ± 0.25 ns in the MOT 
for the 6s level lifetime. Both results are consistent with 
each other and the average gives 45.57 ± 0.17 ns. 

Figure compares the obtained 6s level lifetime with: 
(a) a previous measurement by Marek et a/.,[25| (b) and 
(c) are ab initio many body perturbation theory calcula- 
tion by Safronova et aZ.|26| and W. R. Johnson et aLj^j 
and (d) is a semiempirical prediction using a one electron 
effective potential model by C. E. Theodosiou et al. 28]. 



VI. CONCLUSIONS 

The consistency of lifetime results using two differ- 
ent sources of atoms demonstrates the understanding 



TABLE II: Error budget and corrections for the 6s level life- 
time measurement in the vapor cell 



% error 


Time calibration 


±0.01 


TAC/MCA nonuniformity 


±0.11 


Radiation trapping 


±0.01 


Imperfect laser turnoff 


±0.07 


Magnetic Field 


±0.17 


Background slope 


±0.22 


PMT response 


±0.21 


Bayesian error 


±0.17 


Statistical error 


±0.38 


Total 


±0.56 



TABLE III: Error budget for the 6s level lifetime measure- 
ment in the MOT 









Experiments 






MOT 




present work 1 




vapor cell 








average 




(a) 








(b) ♦ 




Theory 




(c)* 






(d)* 







41 42 43 44 45 46 47 48 49 50 51 
Lifetime [ns] 

FIG. 9: Comparison of the present measurement of the 6s 
level lifetime with previous experimental results (a) by J. 
Marek et al. .|25| and theoretical predictions by (b) M. S. 
Safronova et aL(26| (c) W. R. Johnson et al. |27| and (d) C. 
E. Theodosiou et al. 28] 
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achieved in the systematic effects present in this tech- 
nique. The MOT source advantages as far as density, 
temperature, and radiation trapping compared to a room 
temperature cell; however the intrinsic magnetic gradi- 
ent has been a source of concern. This work sets a lower 
limit on their influence by doing the measurement with 
the two atomic sources. Statistics limit the MOT results, 
while the radiation trapping is the limiting factor in the 
glass cell. This measurement improves previous determi- 
nations by a factor of thirty. 

Our measurement of the lifetime of the first excited 
state in the s manifold of rubidium has reached a pre- 
cision of 0.3 % which permits a very careful tests of the 
calculated matrix elements that contribute to the lifetime 
of the 6s state. These ab initio MBPT calculations repre- 
sent the state of the art in atomic structure calculations 
and require a very thorough understanding of the elec- 
tron correlations and many other subtle effects that add 
and subtract in the expansions. The same techniques ap- 
plied for rubidium are used in cesium and the agreement 
between experiment and theory enhances the confidence 
in the established calculation techniques to provide the 
necessary information such as matrix elements for the 
extraction of weak interaction parameters from the PNC 
measurements in cesium and in future measurement in 
other alkali atoms. 



N 2 „ Ni 
N 2 = + B 2 — 

t 2 n 

N 3 N 
N 3 = l+e 3 N 4 + B 3 ^ 

t 3 n 



eiN 3 , 



(10) 



where Ni represents the number of atoms in state i with 
Nt + N 2 + N 3 + N4 — N the total number of atoms, r, is 
the lifetime of state i, ej stands for the excitation to state 
i, and Bi is the branching ratio from state 1 to state i. 

We take the steady state solution of Eq. ^5] as the 
initial state for the decay. The solution of Eq. ^] with 
ei = e 3 — corresponds to the case of a perfect on/off 
ratio. 

To calculate the perturbation introduced by the imper- 
fect turnoff we use the unperturbed solution for N 3 (t) on 
the other two equations in the following way 



Nt 
N 2 



Ni 

Ti 
N2 
T 2 



f(t)elN 3 



B- 



Ni 



(11) 



with f(t) the second step laser turnoff which is well repre- 
sented by a Gaussian function (At « 5 ns) plus a constant 
level 
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VII. APPENDIX 

We calculate the effect of the imperfect laser turnoff by 
including it in the rate equations as a small perturbation 
and compare the result to the unperturbed solution. For 
simplicity we rename the 6s, 5P 3 / 2 , 5P\/ 2 and 5s levels 
as 1, 2, 3, and 4 respectively. The rate equations with 
the two step excitation on are 

Ni 

Ni = + ei N 3 

Ti 



f(t) = (l-a)exp(-^j +a, (12) 

where a is the on/off ratio. 

We use the analytic solution of Eq. ^2 for N 2 (t) to 
generate a decay sample that we process in the same 
way as real data to obtain a lifetime for the 6s level. 

We include the following parameters in the model: For 
the lifetimes of the 5p levels we use the results 0, El El 
(as in Table P) and for the branching ratios Bi we use 
the theoretical values calculated by Safronova et aLpCjl] 
We take e 3 — 1/t 3 and to obtain ei we compare the 
observed number of counts when the excitation reaches 
steady state to -/V 2 (0) in the model. 

The most likely parameter configuration consistent 
with our experimental setup gives a lifetime correction of 
-0.11 ns. We perform a sensitivity analysis for the model 
parameters to obtain an uncertainty in the lifetime of ± 
0.1 %. 
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